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Abstract. To every Hermitian vector bundle with connection over a compact 
Riemannian manifold M one can associate a corresponding connection Lapla- 
cian acting on the sections of the bundle. We define analogous combinatorial, 
metric dependent Laplacians associated to triangulations of M and prove that 
their spectra converge, as the mesh of the triangulations approaches zero, to 
the spectrum of the connection Laplacian. 



1. Introduction 

Let M denote a compact Riemannian manifold. Combinatorial analogues of the 
Laplacian on M have been considered in various contexts since the appearance of 
[9j. In the 1970s Dodziuk defined combinatorial, metric dependent Laplacians asso- 
ciated to triangulations of M and proved that under certain technical assumptions 
the spectra of these operators converge to the spectrum of the Laplacian acting on 
functions on M as the mesh of the triangulations tend to zero, cf. [6]. This result 
was extended to the Laplacian acting on differential forms in 7J and used in one of 
proofs of the Ray-Singer conjecture, cf. [12 . 

A natural generalization of the Laplacian on M is the connection (also called 
rough) Laplacian in whose definition the exterior d is replaced by the covariant 
derivative V associated to a connection on a Hermitian vector bundle over M . Not 
much seems to be known about the spectra of such operators in general (but see 
[3j). When the bundle is trivial of rank 1 (and usually in a noncompact setting) 
such connection Laplacians are also known and studied in the mathematical physics 
literature, e.g. jiSj, as magnetic Schrodinger operators. 

In this note we propose two definitions of combinatorial analogs of such con- 
nection Laplacians. Both rely on replacing the wedge product of differential forms 
by a graded commutative but nonassociative cup product on simplicial cochains, 
which is probably part of folklore, and whose definition can be found for example 
in In the first discretization scheme, which is applicable only for trivial bun- 
dles, we define combinatorial Laplacians corresponding to a coboundary operator 
twisted by a 1-cochain defined using the above cup product. The convergence of 
eigenvalues in degree zero (i.e. for the Laplacian acting on sections only rather 
than on vector bundle valued forms) in this case is the content of Theorem 13.51 
The second construction of discrete connection Laplacian uses the well known fact, 
cf. |15j . that any bundle with connection can be embedded into a trivial bundle 
in such a way that the connection is induced by the standard connection on the 
trivial bundle. Then the idea is to discretize this "classifying map". We establish 
the corresponding convergence of spectra in Theorem 14.41 
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We mention that our approach leads to a natural generalization of the discrete 
magnetic Laplacians on graphs studied in and [14] to simplicial complexes. 
One can expect that an analogous convergence of spectrum distribution functions 
result for such twisted Laplacians can be proved also for noncompact universal 
covering manifolds, using the method of [5]. Finally, we note that a quite different 
discretization scheme for the degree zero case can be found in [10] and |llj . 



In this note M will always denote a closed oriented Riemannian manifold of di- 
mension N. Let i7(M) be the vector space of all smooth complex valued differential 
forms on M and L^(A(M)), the space of square integrable forms on M with respect 
to the volume element induced by the metric. Let K he a smooth triangulation of 
M and C{K) be the vector space of (oriented) cochains of K with complex coeffi- 
cients. In what follows, ||w|| will always denote the L^-norm of the form uj and ||c||, 
the norm of the cochain c given by the canonical inner product on cochains with 
respect to which simplices are orthogonal. 

Recall that for every q > the corresponding de Rham map : il'^{M) 
(K) is given by R^{u}){a) = J^to, where to G ^'(M) and cr is a g-simplex. When 
q = the de Rham map is defined to be simply the evaluation of a function on the 
0-skeleton of K. 

The Whitney mapping goes in the reverse direction and is defined as follows. 
Let — fivi denote the barycentric coordinate function corresponding to a vertex 
Vi and let a = [vq, ....Vq\ be a g-simplex, where q> Q. We set 



'i=0 

If g = we set [vo) ~ Mo- This expression defines a map: C'^{K) 
L^(A'^(M)). We note that the image of the Whitney map does not consist of 
smooth forms and summarize its basic properties in the following 

Proposition 2.1. i) On the complement of the {N — \)-skeleton of K one has 
W^d^ = dW^ , where d^ denotes the simplicial cohoundary operator of K. 

(ii)Let ja- denote the inclusion map: a ^ M and j* the corresponding pull-back 
map on forms. Then the values ofj*W^c (fore an arbitrary cochain) depend only 
on the value of c at a. Therefore is well-defined on the image ofW^ and one 
has R^W^ = id on C{K). 

(Hi) The support ofW^{a) is contained in the closed star of a. 

For proofs the reader is referred to [20] or [6] . 

The injectivity of the Whitney map allows us to define a new, metric dependent 
inner product on the space of cochains, namely, for ci, C2 S C!{K) we set (ci, C2)w — 
{W^ ci, W^C2)l^. We will refer to this inner product as the Whitney inner product. 



3.1. A commutative cup product. We first let K be any finite simplicial com- 
plex. 

Definition 3.1 ( [4], [1], [21]). For two oriented simplices cri £ Cp{K) and (72 G 
C^{K), set (Ti U (72 = unless ai and (72 meet at precisely one vertex and span a 
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{p + (z)-siniplex t, in which case define cri U (72 = ^{o'l, ^2) (^p^'g^iy T , where the sign 
e(o'i,o'2) = ±1 is determined by the equation orientation{ai).orientation{a2) = 
e(f 1, 02) ■orientation{T) . 

One checks that this defines a graded commutative non-associative bilinear op- 
eration on C{K) with respect to which the coboundary operator is a derivation. 

In the case when K is the underlying complex of a smooth triangulation of a 
compact manifold M it was observed in [5] ( see also [T] and [H] for more details) 
that this cup product admits an alternative description in terms of the Whitney 
and De Rham maps. Namely, for a, 6 S C{K) one has 

aUb = R^iW^aAW^b) (3.1) 

Now let A he a. fixed real smooth 1-form. Consider the twisted exterior differential 
dA — d + iA A . and the corresponding Laplacian A a = d\dA + dAd\. We use the 
above to define discrete analogues of these operators. 

Let a be a fixed 1-cochain. Define the twisted coboundary operator associated to 
it as d^ = d^ +ia\J ., where d^ is the usual coboundary of the simplicial complex 
K. Define the twisted discrete Laplacian by — {d^^* d^ + d^ (d^)*, where 
the adjoint is taken with respect to the inner product given by the Whitney map. 

From now on we consider only triangulations K which are subdivisions of a 
fixed triangulation and whose fullness is bounded away from (see or [7] for a 
definition). Then the cup product on cochains introduced above approximates the 
wedge product on forms according to the following result established in [21 .(The 
case when uj2 — I was proved in [6].) 

Theorem 3.2. Let a be a simplex in K , p - any point in the interior of a and 
Xi,...,xm - local coordinates around a. Let loi,uj2 G 51(M). There exists constant 
C > independent of L0i,u)2,cr and K such that 



W'^ {R'^LOi U R^UJ2) (p) - A W2 {p) I 





duj2 




duji 


^sup \uji \ sup 


dxi 


+ sup W2 sup 


dxi 



(3.2) 



Here the suprema are taken over all partial derivatives and all points in a. 

Let us recall the definition of Sobolev spaces of vector bundle valued forms. Let E 
be a Hermitian complex vector bundle over M and let il{M, E) denote the smooth 
differential forms on M with values in E. We fix a Hermitian connection on E which 
defines a covariant differential V : n*{M,E) n*+^{M,E). Let H''{A{M,E)) be 
the completion of U{M, E) with respect to the norm \\uj\\r = /^^ + AE)^ij-'\'^dvol, 
where A^; = V*V -I- VV*, V* being the formal adjoint of V, is the connection 
Laplacian associated to E and V. 

When E is the trivial fine bundle we write i7''(A(M)) for iJ''(A(M, E)). In this 
case, integrating the estimate in Theorem 13 . 2 1 with 0^2 = 1 and using the standard 
Sobolev inequality (see e.g. [16]) one obtains for every r > + 1 



(3.3) 
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3.2. Convergence of eigenvalues. We set a ~ R^A and use the results above 
to compare the operators dA and . It turns out that the analogue of the identity 
from Proposition [23] (i) holds approximately in the twisted setting. 

Lemma 3.3. For every uj G r2(Af ) and every point p G M we have 

\W^d^R^uj - dAW^R^uj\^ < CA,u>h, (3.4) 

jW^^df i?^cj - W^R^dAUj\j^ < CA,ujh, (3.5) 
where Ca.uj is a constant depending only on A and lu and their first derivatives. 
Proof. In order to show p.4p . we write 

W^^df - dAW^R'^uj = {R'^iA U R'^Lj) ~ iA A W^R^uj 

= W'^ [R'^iA U R'^uS) -iAhuj -iAh {W^R^u - u) 
and then apply Theorem 13.21 and the basic estimate p.3p . Similarly, observe that 

W^d^R'^w - W^R^dAW = {R^iA U R^uj) - W^R^{iA A uj) 

= W'^ {R^iA U R'^lo) -iAAuj + iAhio- R^ {iA A u). 
This easily implies (13. 5p . □ 

Note that dA extends to a map: H^{A{M)) L^(A(M)), denoted also by dA- 
Proposition 3.4. i) For every lo G il(Af ) we have 

WW^d^R'^uj - dAU}\\ < CA^uh, (3.6) 

where Ca,uj is a constant depending only on A and oj and their derivatives, 
ii) For every c G C^{K) we have 

\\W^d'^c-dAW'^c\\ <CA\\W^c\\^h. (3.7) 
As above, Ca depends only on A and its derivatives of order up to r. 

Proof. The estimate p.6p is an easy consequence of p.Sp . Observe that the image 
of C°{K) under the Whitney map consists of continuous piecewise linear functions 
and therefore is in (A(M)). Denoting by for the r-th Sobolev norm on a 
closed A^-simplex a computed using local coordinates in a neighborhood of cr, we 
obtain 



\\W^d^c~dAW^c\\^ <Y^ I \W^d^c-dAW^c\^dvol< 

E ^V^""^' (lli^^cii;:)' < E CAh' iww^cr.f < CaH^w^cWI 

Above we used the local estimate p.4p with oj = W^c, a version of the Sobolev 
inequality which can be applied to each closed A^-simplex (see, e.g. |2i Theorem 
3.9]) and the inequality ^ < C^i/i^. 

□ 
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It is well-known that Aa with domain all smooth forms is an essentially selfad- 
joint elliptic operator on the compact manifold M and hence its closure, denoted 
also by A a, has purely discrete spectrum. Let Aq < Ai < A2 < . . . be the eigenval- 
ues, repeated according to their multiplicities, of acting on functions and let 
^0 — ^1' !i • ■ ■ < ^dimC"(K} denote the eigenvalues of acting on 0-cochains. 
We have the following convergence of eigenvalues result. 

Theorem 3.5. Let j < dimC'^(_ftr). There exist positive constants Ca and Cj^A 
such that following inequalities hold. 



(\[^ + CAh\ 



Proof. We begin by proving the second inequality. It follows from p.7p that for 
every 0-cochain c and some constant Ca > depending only on A and its first 
derivatives one has 



[W'^d^c - dAW^cW < CAh (\\W^c\\ + IUaW^ 



Thus we obtain 



||W^^df c|| > (1 - CAh) \\dAW^c\\ - CAh \\W^c\\ . (3.8) 
Using this, the min-max principle (see e.g. [THj) gives 



, ,^ \\W^d^c 

^ ci...,c,-ieC»{K)<:^0,{c.c^)=0, \\W^c\\ 



fc=l,2,... j-1 

. f {l^CAh)\\dAW'<c\\ \ 

> sup mf " ^ - CAh 

ci,...,c,_ieco(A)c/o,<c,cfc>=o. \ \\W^c\\ I 

k—l,2,...,j — l ^ ' 

> (1 - CAh) sup inf %^ - CAh 

/i,...,/,_i6W[C0(K)] feW[c'>iK}]\{Q}. 11/11 

(/,/fc>=o,fc=i,...,i-i 



> (1 - CAh) y/Xj ~ CAh. 

Hence the second inequality holds. 

To establish the first inequality, denote by Vj the vector space spanned by the 
first j eigenfunctions of A^i. Then it follows from (|3.6p that for each f € Vj \ {0} 
and for some constant Cj,A depending only on j and A we have 



{W''d^R^f,W^d!^f) {dAf,dAf) 



< Cj^Ah 



{W^Ri^f^W^R^f) (/,/) 
Then the proof proceeds in exactly the same fashion as in 6 , Theorem 5.7]. 

□ 

4. The general bundle case 

4.1. The universal connection. Now let i? be a Hermitian complex vector bun- 
dle of rank d over M . We denote the space of forms with values in E defined 
using the inner product on £' by i {h.{M),E) and fix a Hermitian connection on E 
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which defines a covariant differential V : Q*{M,E) Q*+^{M,E). We shall need 
the following theorem on the existence of universal connections due to Narasimhan 
and Ranianan 15J, in the form proved in the appendix of [17j . 

Theorem 4.1. There exist a trivial Hermitian bundle V and an isometric embed- 
ding of bundles i : E ^ V such that V = i* o d o where i* denotes the fiberwise 
adjoint taken with respect to the corresponding inner products. 

Let us recall briefly how the Narasimhan-Ramanan theorem is proved as the 
construction will be needed in the sequel. One first shows that the statement is 
true for trivial bundles. Then one chooses a finite open cover {Ui} of M such 
that E restricted to every Ui is trivial and a partition of unity {tpi} subordinate 
to this cover such that J^i i'f = 1- If *i denotes the already constructed isometric 
embedding from Elu^ to a trivial bundle Vi one defines 



where V = 0j Vi . 

In what follows, we will assume that we have fixed a triangulation Kq fine enough 
so that there exists a finite open cover {Ui} as above with the additional properties 
that (1) each Ui is a subcomplex of Kq, and (2) E restricted to Ui is trivial, and 
we will consider only subdivisions K of Kq. We will also assume that the map i is 
constructed using such an open cover. 

4.2. Twisted de Rham and Whitney maps. We first define the space of cochains 
on which the discrete connection Laplacian will act. We denote the induced map 
from L (A(A/), E) to L (A(Af), V) also by i. We would like to define cochains on 
K with values in the bundle E. To this end, fix a reference point po- in the interior 
of each simplex a in K. 

Definition 4.2. Let C{K,E), the twisted cochains with values in E, denote the set 
of all maps from the set of all simplices in K to the total space of E which assign 
to each simplex a a vector in the fibre E.p^ . 

This is a vector space with an inner product induced by the Hermitian structure 
on E. The restriction of i to fibers defines a map, denoted by i^, from C{K, E) to 
C{K,V) which we identify with the cochains taking values in C", where n is the 
rank of the trivial bundle V . 

We shall also need the following notation. Let : C(Ui, C^) C{Ui, C") be 
the operator given by {I^c)p = Y^^R^ {ips) U c^. Here ips are the entries of the 
matrix of i and Cp the components of the vector valued cochain c. We define the 
operator ^-f by ^-f Cp = {il^i) U Cp. Finally we define : C{K, E) C{K, V) 
by I^c — J2i ^i^^fc. Now define an operator acting on C{K, E), the discretization 



In order to compare V and V we have to introduce the appropriate twisted 
Whitney and de Rham maps. We set = i*W'^I^ and = (/^)* Here 
and are the usual Whitney and de Rham maps but now acting between the 
spaces of vector valued cochains and differential forms, C{K, V) and L (A(Af ), V), 
respectively. 

It is easily seen that is injective on 0-cochains. More generally, we have: 




(4.1) 
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Lemma 4.3. For sufficiently fine triangulations K the map is injective on 
C^iK^E) for allq. 

The proof will be given in the next section. We can now define the corresponding 
Whitney inner product on C{K,E) as in Section 2 and the discrete connection 
Laplacian = (^^)* + (^^)*' where the adjoint is taken with respect 
to this inner product. Let Aq < Ai < A2 < . . . denote the eigenvalues repeated 
according to their multiplicities, of (the closure of) the connection Laplacian 
Ae = V*V + VV* acting on sections of E, and let A^ < Af < . . . < \fi^c"{K e) 
denote the eigenvalues of acting on 0-cochains with values in E. We have the 
following convergence of spectra result. 

Theorem 4.4. There exist positive constants and such that following in- 
equalities hold for j < dim C° {K, E) . 

(.fxf + C'hY 
Af - C)h < A, < ; ^.^ ^ . 

5. Proof of Theorem 14.41 

5.1. Estimate from belov^r. In this subsection we estimate from below the eigen- 
values of the connection Laplacian in terms of the eigenvalues of its discretizations, 
i.e., we prove the first inequality in Theorem 14. 41 We denote the projection ii* by P 
and introduce the "almost" projection operators — ij^) ■ next lemma 
states that the Whitney and de Rham maps approximately interwine P and P^ . 

Lemma 5.1. For r > ^ + 2, every form uj E i}{AI,V) and every vector valued 
function f G ^l'^{M, V) one has 

\\W^P^R^uj - W^R^Puj\\ < ah\\w\\r, (5.1) 

WW'^P^R^UJ - PW^R^uj\\ < ah\\w\\r, (5.2) 

P'^R'^f = R'-^Pf. (5.3) 
Here Ci denotes a constant depending only on the map i and its derivatives. 

Proof. We first show that (|5.ip holds. In the course of the proof, various constants 
depending on i will all be denoted by Ci. 

We first assume that the bundle E is trivial and identify C{K, E) with C{K, C^) 
and n{M,E) with n{M,C'^)- Let cr be a g-simplex in K and let uj e n^iM^E). 
We write R^ (respectively W^) for both de Rham (Whitney) maps defined on 
CiK,^^) and C{K,V) (respectively on 17(M,C'') and n{M,V)) and observe that 
in view of [71 Lemma 7.22] for every c € C(K, V) one has 

Ci/l^-2g||^||2 < \\W^c\\^ < C2/l^-'«||cf , (5.4) 

where the constants Ci and C2 do not depend on K. From the definitions we have 
I^R^uj - R^iuj = (R'^iips) U R'^LOs - R^iipstOs)). Thus, combining 
(13.21) and (15.41) we conclude that 



< J2 dl^"^ (R'^iips) U R'^iUs - tpsU^s) II + ||W^^i?^ 
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Now (|5.4p implies that 

Wl'^R^'u^R'^iuof <C,e^-''+^\\u\\l. (5.5) 
It is easy to check that {{I^)*c)p — J2s -^^(*ps)'^s- Hence we also have 

||(/^')*i?^CJ - R^i*Uj\\' < C,/l29-^+2||c^||2. (5.6) 

Next, observe that the operators are uniformly bounded with respect to K. 
Indeed, from the combinatorial definition of the cup product and the mean value 
theorem one easily obtains 

Wi'^R'^oj - /^i?^c^f < a/i29-^+2||c^||2 (5.7) 
and is certainly uniformly bounded. We use this fact in the following estimate. 

< ||/^(/^)*i?^'c^ - /^P^rc^f + ||/^i?^z*u; - i?^u*wf < ah^'^-^'+^Ml. 

Applying (|5.4p again, we see that 

WW^pKj^K^ - W^R^Puj\\ < ah\\uj\\r. 

Now let E be an arbitrary bundle. We use the open cover {[//} and the partition 
of unity {-0;} with the properties described after Theorem 14. II and define Pi = iiii- 
Then, using (|4.ip . we find that P = Yli k'^i'^i'^ki*k — J2i''-i^f'''i- the other 
hand, P^ is not equal to J2i ^i^ i^f)'^i^i^)* due to the non-associativity of our 
cup product. However, according to Theorem 13.21 this product is approximately 
associative, therefore we can replace P^ by J2i ^i^ ™ ^^'^ following 
estimate. 

IIP^P^c. - P^Pc^ll < J2 )'(/f - R^'ii^Pfi^uW + 0{h) 

I 

< E f - i?^^r-)|| + )^p%*c. - 

Above we used (j5.5p and (|5.6p already proved over each Ui, and the inequality 
Wi^ffR^'co - R'^^fi^f < C,h^'^-''+'\\iu\\l 



which is proved in exactly the same manner as ()5.5|) . Thus ()5.ip is shown to hold 
for arbitrary E. 

The estimate (|5.2p can be easily derived from (|5.ip . Finally, the equality (|5.3p 
follows directly from the definitions. □ 

We use this lemma to deduce that the combinatorial connection approximates 
the smooth one in the appropriate sense. 

Proposition 5.2. For r > + 1, every section f G f2''(Af, E) and every 0-cochain 
c e C{K, E) one has 

||M?^V^P^/-V/|| <Cl/i||/||.. (5.8) 
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Proof. We estimate as follows, using (|5.3p . (|5.2p and (13. 3p . 

||W?^V^^^/ - V/ll = ||rW^P^rf^P^i?^i/ - t*dif\\ 

^ \\t*W^P^R^dif - i*dif\\ < IKVK^P^i?^ - P^M^^i^^)d^/|| 

□ 

The first inequality in Theorem 14.41 follows as in the proof of Theorem 13.51 

5.2. Estimate from above. In this subsection we estimate from above the eigen- 
values of the connection Laplacian in terms of the eigenvalues of its discretizations, 
i.e. we prove the second inequality in Theorem 14.41 This will follow from Proposi- 
tion [5l4] below. We will need the following strengthened version of Lemma [5Tl 

Lemma 5.3. For r > + 2 and every form uj € ^1{M, V) one has 

\\dW^P^R^uj - dW^R^Puj\\ < Cih\\w\\r, (5.9) 

\\dW^P^R^uj - dPW^R^uj\\ < ah\\w\\r. (5.10) 
Above Ci is a constant depending only on the map i and its derivatives. 

Proof. The argument is similar to the proof of Lemma |5. II but now one has to use 
in addition the derivation property of d^ . The details are omitted. □ 

We now use this lemma to derive the following proposition which is the analog 
of part ii) of Proposition 13.41 in this situation. 

Proposition 5.4. For every 0-cochain c G C{K,E) one has 

WW^V^c - VW^^cll < C'h{\\W^c\\ + ||VVK^c||). (5.11) 

Proof. We first prove the analogous estimate in the case when c = R^ f for some 
smooth section /. Using (|5.2p . (|5.10p and (|5.3p one has for r large enough: 

\\W^V^R^f - VW^R^fW = WeW^P'^R'^dif - i*dPW^P^R'^if\\ < 
WiW^P'^R^ _ PW'^R^)dif\\ + WdW^P^R^if - dPW^R^tfW < C'^hWfWr. 

Now we note that (|5.3p implies that c — R^W^ and proceed exactly as in the proof 
of part ii) of Proposition 13.41 The only difference now is that the image of 
docs not consist of piecewise linear sections of E anymore. However one can pass 
from the r-th to the first Sobolev norm on each TV-simplex simply by a repeated 
use of the Leibniz rule. □ 

It remains to prove Lemma 14.31 Observe that according to (|5.7p it suffices to 
prove that i*W^i^ is injective. Suppose that cq G C'^{K,E) is a nonzero cochain 
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supported on a single g-simplex a with reference point po-. A direct computation 
in local coordinates shows that i*W^ cq is not identically zero. Now let c be 
arbitrary nonzero cochain and suppose that i*W^i^c = 0. Suppose also that c 
is not zero at a simplex a. Then Proposition 12.11 will imply that j*{i*W^i^c) = 
i*j*{W^i^c) is not zero, a contradiction. 
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